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Abstract 

Cooperative hybrid-ARQ (HARQ) protocols, which can exploit the spatial and temporal diver- 
sities, have been widely studied. The efficiency of cooperative HARQ protocols is higher than that 
of cooperative protocols, because retransmissions are only performed when necessary. We classify 
cooperative HARQ protocols as three decode-and-forward based HARQ (DF-HARQ) protocols and two 
amplified-and-forward based (AF-HARQ) protocols. To compare these protocols and obtain the optimum 
parameters, two unified frameworks are developed for protocol analysis. Using the frameworks, we can 
evaluate and compare the maximum throughput and outage probabilities according to the SNR, the relay 
location, and the delay constraint for the protocols. 
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I. Introduction 

In wireless mobile communication systems, various diversity techniques, such as time, fre- 
quency and spatial diversity, have been investigated to achieve spectrally efficient and reliable 
communications over fading channels [HI -[El. Among them, cooperation diversity techniques, 
which have been widely studied at present, can provide spatial diversity by cooperating between 
users. In cooperative communication systems, distributed antennas of different users are formed 
as a "virtual array" by sharing their antennas and time/frequency resource to achieve the spatial 
diversity. To efficiently obtain the diversity, there have been many approaches for developing 
cooperative protocols, especially for half duplex terminals since a full-duplex terminal is of 
less interest to practical applications [[6l-[[9|[. In [[3, the selection relaying protocol based on 
channel measurements between terminals and the incremental relaying protocol based on a 
limited feedback from the destination terminal were proposed. These protocols were designed for 
time-division multiple-access (TDMA) systems in which either the source or the relay terminal 
transmits signals at a time slot. In [8J, three cooperative protocols for amplified-and-forward (AF) 
and decode-and-forward (DF) modes were proposed, in which the source and the relay terminals 
transmit signals simultaneously to improve performance. However, these protocols are performed 
poorly in high spectral efficiency scenarios, because twice the bandwith of direct transmission 
is required for a given rate. Thus, in [9J, the dynamic decode and forward (DDF) protocol was 
proposed, in which the relay listens to the source until it can successfully decode the information 
bits before transmitting. The DDF protocol achieves the optimal diversity-multiplexing (D-M) 
tradeoff when the multiplexing gain, re, is in the region of < rg < 0.5. 

Automatic repeat request (ARQ) is a common technique used to make a wireless link reliable. 
The cooperative protocols can adopt the ARQ technique by exploiting feedbacks from the relay 
and destination terminals. Since the relay or source retransmits only when the destination terminal 
wants to, the efficiency of cooperative ARQ protocols can be improved. This performance 
improvement of cooperative ARQ protocols has been shown in literature, such as in [[3, [[T0|. 
and fn}. In [Q, it was shown that the incremental relaying protocol, which can be viewed as an 
extension of a hybrid- ARQ into a cooperative context, outperforms the fixed relaying protocol 
in terms of its outage behavior. In [llOil . a DDF based ARQ protocol was proposed for two 
cooperating single- antenna terminals and a double-antenna destination terminal. It was shown 
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that it can achieve the optimal D-M tradeoff when the number of retransmissions goes to infinity. 
In [fm . it was shown that the DDF based ARQ protocol can achieve the optimal D-M tradeoff 
in a single user relay channel when the maximum allowable number of transmissions is greater 
than 2. The performance of a cooperative ARQ protocol and that of a non-cooperative ARQ 
protocol were also compared in [11211 . It was shown that a practical cooperative ARQ protocol 
using a convolutional code is better than a non-cooperative ARQ protocol. 

Inspired from these performance improvements, extended versions of previous cooperative 
protocols combined with ARQ or hybrid- ARQ (HARQ) scheme were proposed and analyzed 
in lfT3l - [[T5l . In [fT3l . three cooperative ARQ protocols, which combine the incremental relaying 
with the selection relaying, were proposed. The outage behaviors of the three protocols were 
shown for a simple ARQ scheme without packet combining. In [1141 and [|15i[ . cooperative HARQ 
protocols for multiple relays were proposed for various HARQ techniques. The upper bound of 
an incremental redundancy based protocol was developed in [fTSl . The achievable throughput 
of these cooperative ARQ protocols changes according to channel and environment such as the 
relay location, the path loss, and signal-to-noise ratio (SNR). Thus, it is difficult to compare these 
cooperative ARQ protocols using previous analysis methods such as the outage behavior and 
the D-M tradeoff. In the D-M tradeoff analysis, the relay location, which is one of the dominant 
factors in determining performance, is not considered. Some protocol parameters, such as the 
initial transmission rate and the maximum number of transmissions (or the delay constraint), 
cannot be optimized because the D-M tradeoff provides a fundamental but only asymptotic 
performance. Outage behavior is also highly affected by the initial transmission rate. Thus it is 
necessary to develop unified frameworks that can analyze and compare protocols considering 
each protocol's characteristic. 

In this paper, two unified frameworks are developed. One is for DF based cooperative HARQ 
(DF-HARQ) protocols and the other is for AF based cooperative HARQ (AF-HARQ) protocols. 
In this paper, the HARQ technique mainly considered is the incremental redundancy HARQ 
(IR-HARQ) technique due to its better performance over the Chase combining HARQ scheme 
[fT6l . In the DF-HARQ protocols, only IR-HARQ scheme is used. On the other hand, in the AF- 
HARQ protocols, both the IR-HARQ and the Chase combining schemes are considered because 
the relay terminal just forwards the amplified packet to the destination terminal. In this paper, we 
consider the extension of cooperative protocols or the cooperative ARQ protocols into cooperative 
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HARQ protocols. These were then classified into three DF-HARQ protocols and two AF-HARQ 
protocols. In the DF-HARQ protocols, Protocol 1 includes the selection and incremental relaying 
DF protocol in Q, and protocol 1 in [fT3l . Protocol 2 is an extended version of protocol 2 in 
dHl. Protocol 3 includes the protocol with a single relay in |fT5l and the DDF protocol when 
Te ^ 0.5 in [9] and UllL In the AF-HARQ protocols, protocol 1 includes the selection and 
incremental relaying AF protocol in fj], and protocol 3 in |[8l. Protocol 2 includes protocol 3 in 
jH and the NAF protocol when re ^ 0.5 in Q- Furthermore, the performance of each protocol 
is evaluated under two different power constraint scenarios. One is peak power constraint, where 
each terminal uses the same power for a transmission. The other is step power constraint where 
the transmission power for each step is preserved. Thus, the analytical framework provided in 
this paper is quite meaningful because: i) it provides a unified approach for the analysis of 
cooperative HARQ protocols; ii) it provides actual performance comparison among different 
protocols, which can be only done in an asymptotical manner with previous approaches; and iii) 
such an analysis, considering the characteristic of each protocol, can be used to optimize the 
performance of each protocol according to channel and environment. 

The rest of this paper is organized as follows. In Section UIl the signal model and the 
cooperative HARQ protocols are described. In Section [nil the analytical frameworks are de- 
veloped for the DF-HARQ and the AF-HARQ protocols. In Section |IVl the outage probabilities 
of terminals participating in each protocol are analyzed using the proposed framework. The 
maximum achievable throughput of the cooperative HARQ protocols are obtained and compared 
in Section |Vl Finally, the concluding remark is given in Section |Vll 

II. System model and protocol description 

A. Protocol Description 

In this paper, a cooperative communication model with three half-duplex terminals is con- 
sidered. It consists of a source terminal, a destination terminal, and a relay terminal. The 
source terminal transmits data to the destination terminal with the aid of the relay terminal. 
The destination and the relay terminals transmit ACK/NACK signals for retransmissions. In this 
subsection, we describe three DF-HARQ protocols and two AF-HARQ protocols. The DF-HARQ 
protocols are comprised of two steps and can be described as follows: 
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• Protocol 1: In Step 1, the source terminal broadcasts a packet to the destination and the relay 
terminals until either one of the terminals successfully decodes the packet. If the destination 
terminal successfully decodes the packet, the source terminal broadcasts the next packet at 
the next time slot. If the destination terminal fails but the relay terminal decodes the packet, 
Step 2 begins. In Step 2, the relay terminal transmits the packet to the destination terminal 
until the packet is successfully decoded. 

• Protocol 2: In Step 1, the source terminal transmits a packet to the relay terminal only. 
If the relay terminal successfully decodes the packet. Step 2 begins. In Step 2, both the 
source and the relay terminals transmit the packet using a space time code to the destination 
terminal until it successfully decodes the packet. 

• Protocol 3: The procedure of Step 1 is equal to that of protocol 1 and the procedure of 
Step 2 is equal to that of protocol 2. 

The AF-HARQ protocols are comprised of two steps and can be described as follows: 

• Protocol 1: In Step 1, the source terminal broadcasts a packet to the destination and 
relay terminals. If the destination terminal decodes the packet successfully, then the source 
terminal transmits the next packet. If the destination terminal fails, the relay terminal 
amplifies the signal and transmits it in Step 2. After that, if the destination terminal fails 
to decode the packet, the next round starts with the same packet. 

• Protocol 2: Step 1 of protocol 2 is equal to that of protocol 1. In Step 2, both the source 
and the relay terminals transmit packets to the destination terminal using a space time code. 

Figures [T] and [2] respectively depict the DF-HARQ and the AF-HARQ protocols. In these fig- 
ures, 'S^ D,R' represents that the source terminal (S) is broadcasting a packet to the destination 
terminal (D) and the relay terminal (R). 'R^D' represents terminal R transmitting to terminal 
D. 'S,R^D' represents both terminal S and terminal R transmiting to terminal D using a space 
time code. Also, 'x(o)' and 'x(x)', respectively, denote that the terminal has achieved decoding 
successfully and has failed decoding, where x ^ 

B. Signal Model 

It is assumed that a packet is composed of one codeword and the channel state information 
is perfectly known at the receiver. Also, the ACK/NACK signaling is assumed to be error-free 
for simple analysis. The resource spent for signaling ACK/NACK is also ignored since it is 
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typically quite small compared to that used for data transmission. Let Xa [m] denote the m-th 
transmitted packet from the terminal a G {S,D,R]. Then, the m-th received packet in terminal 

13 = {S, D,R}\a can be written as 



where Eg is the symbol energy and /i^^^ [m] denotes the channel coefficient between a and (3. 
Here, the channel is assumed to be a block fading channel. It has independent and identically 
distributed (i.i.d.) circularly symmetric complex Gaussian random variable with distribution 
CA/'(0, 1). Also, n/3[m] is an i.i.d. circular complex additive white Gaussian noise (AWGN) 
at terminal /3 with CJ\f{0, Nq). At the source terminal, the information of b bits is encoded using 
the channel code with codebook C E C^^^ of length LM over the complex numbers. The overall 
codeword is divided into M blocks of length L symbols. The i-th codeword C, obtained from 
the codebook C is transmitted at each transmission. Each packet is transmitted during r (sec) 
over bandwidth W (Hz) and it is assumed that a packet has one codeword so that a codeword 
has L ^ Wt degrees of freedom. The initial transmission rate (spectral efficiency) is r = b/L 
(bit/sec/Hz) and the initial transmission rate is selected based on the long-term average signal 
to noise ratio (SNR). The maximum number of transmissions is M, i.e. the maximum number 
of retransmissions is M — 1. 

In the DF-HARQ protocols, if the relay terminal decodes the packet successfully, the relay 
transmits a different codeword than that of the source terminal as follows. The information of b 
bits is encoded using the channel code with codebook C G C" where n = 2LM. The codeword 
set Cs = {Ci, ■ ■ ■ , Cm} is used at the source terminal for transmission and the codeword 
set Cr = {Cm+1, ■ ■ ■ ) C2j\/} is used at the relay terminal. In the AF-HARQ protocols, the 
source terminal transmits different codewords at every transmission. However, the relay terminal 
forwards the amplified packet to the destination terminal. Thus, the amplified packet received 
from the relay terminal and the packet received from the source terminal are combined using 
the Chase combining method at the destination terminal. 

III. Analytical frameworks for cooperative hybrid-ARQ protocols 

In [[TtI and [fTSl, the throughput of a point-to-point link using an HARQ scheme is analyzed. 
The event that a terminal stops transmitting the current codeword is recognized to be a recurrent 
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event [fT9ll . A random reward TZ is associated to the occurrence of the recurrent event. If successful 
decoding is achieved within M transmissions, IZ = r b/s/Hz. If successful decoding is not 
achieved, 7^ = b/s/Hz. Using the renewal reward theorem llT9l , the system throughput is 
obtained as 

where T is the number of blocks spent to transmit a single packet. This equation can be directly 
used for cooperative HARQ protocols. Let Pout be the outage probability of cooperative HARQ 
protocols, then (|2l) is given as 

= (3) 

^ E[T] 

Thus the throughput of cooperative HARQ protocols can be obtained by computing po^t and 
E [T]. In this section, two unified frameworks for the DF-HARQ and the AF-HARQ protocols 
are developed to obtain pout and E [T]. 

A. A framework for DF-HARQ protocols 

In this subsection, a unified framework for the DF-HARQ protocols is developed using state 
transition diagram approach. Let 5*^ {Sm) be the event denoting successful decoding (decoding 
failure) at the m-th transmission and q (m) be the probability of having successful decoding at 
the m-th transmission, defined as 

q{m) = ¥1 {Si,S2,- ■ ■ ,Sm-i,Sm) ■ (4) 

Then, 

= Pr (5'i,5'2,- ■■ ,5"^) (5) 
is the probability of decoding failure with m received packets and q (m) can be rewritten as 

g(m) = p{m — 1) — p (m) . (6) 

The probability functions of the destination and the relay terminals in Steps 1 and 2 are defined 
as follows. Respectively, pi (m) and p2 (m) are the error probabilities of the destination terminal 
and the relay terminal in Step 1 when the m-th packet is received from the source terminal. 
The outage probability with m received packets of the destination terminal in Step 2, while the 
relay decoded it successfully with n transmitted packets in Step 1, is defined as p3 „(m). Also, 



December 16, 2008 



DRAFT 



JOURNAL OF MbX CLASS FILES, VOL. 1, NO. 13, MAY 2008 

TABLE I 

Transition probabilities for the DF-HARQ protocols 



8 





Transition probabilities 


dm 


pi (m|m ~ 1) P2 (m\m — 1) T 


bm 


pi (m\m — 1) (1 — p2 (m|m — 1)) T 




(1 - pi (m|m - 1)) TV, (m < M) 


Cm 


(1 - pi {M\M - 1)) TV + pi [M\M - 1) p2 {M\M ~1)T 




P3,n (m m — 1) T 




(1 - P3,n (m|m - 1)) TU\ (m < M) 


CM.n 


T 



gi (m) , g2 and g3,„ (m), respectively, denote the probabilities of successful decoding with 
m received packets at the destination and the relay terminals in Step 1, and at the destination 
terminal in Step 2. 

Figure [3] shows the state transition diagram of the DF-HARQ protocols. Here, all branches are 
labeled with the corresponding transition probabilities, multiplied by dummy variables T and U. 
The exponents of T and U, respectively, denote the number of transmission and the throughput 
during the transition. The states are labeled as: starting of one round (Xg); the state when the 
source terminal transmits a packet at the m-th transmission in Step 1 (Am)', the state when the 
relay terminal transmits the packet at the n-th transmission in Step 2 after the relay successfully 
decodes the received packet at the m-th transmission and the end of the round when the 

destination terminal decodes the packet successfully or the number of transmissions has reached 
M (Xe). The transition probabilities for each state are listed in Table HI Here, p{m\m — 1) 
denotes the conditional probability of decoding failure at the m-th transmission, given that the 
receiver has not decoded it until the (m — l)-th transmissions. As derived in Appendix lAl the 
average number of transmissions can be calculated using the state transition diagram as 

M-l M 

E[T] = 1 + Y,P^ ('^) P2 (m) + Y,Pi (^) 92 (n) E [T'] , (7) 

m=l n=l 

where E \T'] = 1 + Ylm=iP'i,n {m). Also, the outage probability of a DF-HARQ protocol is 
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TABLE II 

Transition probabilities for the AF-HARQ protocols 





Transition probabilities 




pi (m|m — 1) r 


bm 


P2 (m m — 1) r 




(1 -Pi (m|m- 1))TU'' 




(1 - P2 (m|m - 1)) ^f/^ m < M 




T 



given by 

M 

Pout = Pi Pi {n) q2 (n) P3,n {M) . (8) 

n=l 

By substituting (|7]) and ([8]) in ([3]), the throughput of the DF-HARQ protocols can be obtained. 

B. A framework for AF-HARQ protocols 

To obtain the throughput and the outage probability of the AF-HARQ protocols, a unified 
framework for the AF-HARQ protocols is developed using a state transition diagram. The error 
probability at the destination terminal is defined as follows. Respectively, pi (m) and p2 (m) 
are the error probabilities at the destination terminal in Steps 1 and 2 when the m-th packet is 
received from the source terminal. Also, qi (m) and q2 (m), respectively, denote the probabilities 
of successful decoding at the destination terminal in Steps 1 and 2. 

Figure m shows the state transition diagram for the AF-HARQ protocols. In the state transition 
diagram, all branches are multiplied by dummy variables T and U as in the case of the DF- 
HARQ protocol. The states are labeled as: starting of one round (X^), the state when the source 
terminal broadcasts a packet in Step 1 at the m-th transmission (Am)', the state when the relay 
terminal relays an amplified packet to the destination terminal in Step 2 at the m-th transmission 
(Rm)', and the end of the round when the destination terminal decodes successfully or the number 
of transmissions has reached M (Xg). The transition probabilities for each state are listed in 
Table ini As derived in Appendix El the average number of transmissions can be calculated using 
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the state transition diagram as 



M 



E\T]= ^ (2m — (m) p2 {m — 1) + 2mpi (m) q2 (m). 



(9) 



m=l 



Also, the outage probability of an AF-HARQ protocol is given by 



Pout=Pl{M)p2 (M). 



(10) 



By substituting Q and (flOl) in (O, the throughput of the AF-HARQ protocols can be obtained. 

Because the two frameworks consist of the error probability of each terminal, the throughput 
and the outage probability of any cooperative HARQ protocol can be obtained from knowing 
of the error probability corresponding to a given protocol at each terminal. Thus, the proposed 
analysis can be used as framework for the analysis of any cooperative HARQ protocol. Also, 
the upper limit of protocol performance can be obtained from an information theoretic measure, 
as well as, the actual performance from the measured or simulated error probabilities using 
practical modulation and coding schemes. 



In this section, the outage probability at each terminal is calculated to obtain the performance 
of the cooperative HARQ protocols using the proposed frameworks in Section Hill Throughout 
our analysis, Gaussian codebooks are assumed to obtain the upper limit of performance. The 
probability functions are approximated using the central limit theorem to make the analysis 
tractable. Later, simulation results will show that the approximation is quite accurate. 

A. Outage probability 

In ifTTl . it was shown that the probability of error can be arbitrarily small if > r, or almost 
one if Im ^ r when a length of the codeword L is large enough. Thus we can assume that 

Pr \error\Irr, > r| = 0, 

1 I m J , ^^^^ 

Pr {error\Im ^ r} = 1, 

and the probability of undetected error is 0. Thus the error probability p{m) of each terminal 
can be replaced by the outage probability of each terminal. The outage probability at the m-th 
transmission is defined as 



IV. Outage probability at each terminal 



p (m)=Pr [S^) = Pr (5*1, 5*2, , S^) 



(12) 



Pr {Im < r) , 



(13) 
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where r is the initial transmission rate and is the accumulated mutual information at the m-th 
transmission. The accumulated mutual information for HARQ systems can be calculated as ^lOl 



m 



/n. = J]J J]7.,s' , (14) 
s=l \s'=l J 

where 7., s' = ^ denotes the average SNR. Here, m' denotes the number of packets 

for the Chase combining, and m denote the number of packets for the code combining in the 
IR-HARQ scheme. For the Gaussian inputs, the mutual information is explicitly given by 

J (7„) = log2 (1 + 7„) . (15) 



B. DF-HARQ protocols 

In this subsection, the outage probabilities pi(m), P2{rn) and Ps^nif^) are calculated for each 
protocol. In the DF-HARQ protocols, the source and the relay terminals transmit different 
codewords at every transmission. 

1) Protocol 1: From (fT2l) . pi{m) and P2{m) are respectively given as 



Pi [m) 



,s=l 



P2(m) = Pr|^5^j(7f)<rj, (17) 

where 7s = ^ |/is,D[s]|^ and 7^ = ^ l^s,/?!^]!^. Because the directly received and the relayed 
packets have different redundancy information, P3,„(m) is given as 

(n m \ 

s'=l s=l J 

where 7^ = ^ |^. To make the analysis tractable, a Gaussian approximation ifTSl can be 

used when m and n axe not small. Then, pa (m), a = 1,2, can be approximated as 

Pa{m)^l-Q ('-^^] , (19) 

V V^^a J 

where /i^ and are the mean and the variance of J (7m)- Also, Ps^ni'm) can be approximated 
as 

ft,„W^i-Qf''7'"f-":M. (20) 
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where /xa and (j| are, respectively, the mean and the variance of J (7^ ) and /X3 and a'^ are, 
respectively, the mean and the variance of J {'y])- Thus the outage probabilities can be obtained 
from the means and the variances of J (7}), J {'Ji), and J {'jf). 

2) Protocol 2: The destination terminal does not receive packets from the source terminal in 
Step 1. Thus pi{m) can be set to 1 and P3,„(m) can be simplified as Ps^m), (i.e., J (7]/) = 0). 
The outage probability P2{m) in protocol 2 is equal to the that of protocol 1. Since both the 
source and the relay terminals transmit packets using a space time code, P3{m) can be given as 

P3 (m) = ^'{j2-J hi + Is) < ■ (21) 

The Gaussian approximated outage probabilities of P2{rn) and P'i{m) can be obtained using (fT9l) . 

3) Protocol 3: In protocol 3, Step 1 is equal to that of protocol 1. Thus pi{m) and P2im) 
of protocol 3 are equal to those of protocol 1. Step 2 of protocol 3 is equal to that of protocol 
2. However the destination terminal of protocol 3 has the n directly received packets from the 
source terminal in Step 1. Thus, p-i,n{m) is given by 

(n m \ 

E J hi) hi +^s)<r]. (22) 

s'=l s=l J 

The Gaussian approximations of pi{m) and P2{m) can be obtained using (fT9l) . Also, the Gaussian 
approximations of ps,n{ni) can be obtained using (|20|) . 

The means and the variances of J(7i), a = 1,2,3, can be calculated as follows. The mean, 
Ha, is obtained as 

A'a=£[J(7ni = / log2(l +7r)^exp(-^)d7r (23) 
--Pl-i'U,: f _ M , (24) 



ln(2) V 2(t2^ 

where Ei is the exponential integral [22J and is the channel variance. It is difficult to calculate 
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the variance of J(7i) exactly. However, if SNR is high, it can be approximated as 



= E [(log,(l + 7r))V/^a 

[(log,(7?))V/^a 



(25) 



(26) 



where C is the Euler number ||22]| . If SNR is very small, it can be approximated as 



(27) 



where e is the base of natural logarithm. For moderate range of SNR, the variance is obtained 
by using numerical integration. 

C. AF-HARQ protocols 

In this subsection, the outage probability of each terminal is calculated according to the AF- 
HARQ protocols. In the AF-HARQ protocols, the relay amplifies a received signal with an 
amplifier gain a^, defined as 



where Eg denotes the average transmission energy at the source terminal and Er denotes the 
average transmission energy at the relay terminal. At the destination terminal, the packets received 
from the source and the relay terminals are combined using the Chase combining method. 
1) Protocol 1: Respectively, pi{m) and P2{rn) are given as 




(28) 




(29) 



and 




(30) 
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where 



Is + 7f + 1 ' 

For the AF-HARQ protocols, also the Gaussian approximation [fTSll can be used when m and n 
are not small. The Gaussian approximated pi{m) is given as 

/ N . ^ I r — (m — l)ui — u\\ 
pAm)^l-Q\ — / ^ =J- , (31) 

where ui and vf are, respectively, the mean and the variance of J (7J + Fi) and u[ and f f are, 
respectively, the mean and the variance of ^(7^). Also, the Gaussian approximated P2{rn) is 
given as 

p.(™)«l-Qf^V (32) 

where U2 and are, respectively, the mean and the variance of J (7J + Fi). 

2) Protocol 2: The source and the relay terminals transmit space time coded signals to the 
destination terminal in Step 2. Thus, F^ can be written as 

= + 2 t'^h 1 • (33) 
Ts + 7f + 1 

By substituting (l33l) into (|29l ) and (|30l) . the outage probabilities of protocol 2 can be obtained. 

To obtain the mean and the variance of the J (71 + Fi), the probability density function (pdf) 
of Fi is required. However, it is difficult to obtain the exact pdf of Fi. Instead, half the harmonic 
mean of 7^ and 7^ is used as in [|2T1l . which is given by 

Ts + Ts 

and the pdf of H is given as 11211 

■ 27 



/«(7)=^e--(^/^)i^o ( 



+ ^e-(-/-)ir,f^y (35) 
V^/V Vv^/ 

where U = -^[7^] + E[y^], V = E['y'^]E['y^], and Ky,{-) denotes the modified Bessel function of 
second kind with the order of w. Respectively, the mean and the variance of J (7^ + Fi), n and 
t>^, can be obtained as 

Ua = E[J {-fl + Fi)] =E[J{X + H)] (36) 
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and 

vl=v[J (7i+ri)] 

=E[J{X + nf]-E[J{X + n)f (37) 

where X = 7I for protocol 1 and X = 2'jI for protocol 2. It is difficult to integrate (l36l) and 
(l37l) into closed forms and we resort to a numerical integration method to calculate (|36l) and 

V. Numerical results 

In this section, simulation is performed and compared with the analysis for the DF-HARQ 
and the AF-HARQ protocols. The cooperative system is assumed to be a one-dimensional linear 
relay network for simplicity. Let the distance between the source and the destination terminals 
be 1 and the distance between the source and the relay terminals be d. Channel is assumed 
to be an i.i.d Rayleigh block fading. The long-term average channel gain of each link is set 
to be gs,D = 1 for the source-to-destination link, gs,R = d^"" for the source-to-relay link, and 
9r,d = (1 ^ for the relay-to-destination link where the path-loss exponent a is set to 4. 
The optimum initial transmission rate is stepwisely searched over ~ 12M with the step size 
of 0.2M. For performance comparison, we consider the peak power constraint, where the source 
and the relay terminals transmit packets with the same power, and the step power constraint, 
where the total transmission power for each step is constrained and the transmission power in 
Step 2 is equally divided to the source and the relay terminals, (i.e., the source and the relay 
terminals transmit packets in Step 2 with half the transmission power of the source terminal in 
Step 1). 

A. Peak power constraint 

In this subsection, we obtain the maximum throughput of the DF-HARQ and the AF-HARQ 
protocols under the peak power constraint. In these figures, 'DFl', 'DF2', and 'DF3', respectively, 
denote the analysis of protocols 1, 2, and 3 of the DF-HARQ protocols, and 'AFl' and 'AF 
2', respectively, denote the analysis of protocols 1 and 2 of the AF-HARQ protocols. Also, 
Es/Nq denotes the long-term averaged SNR of the source signal at the destination terminal. 
Respectively, 'DF simuL' and 'AF simuL' denote the throughput obtained from the simulation 
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of the DF-HARQ and the AF-HARQ protocols. From Figures [5] - |7l it is shown that the analysis 
for the DF-HARQ protocols matches with the simulation results quite accurately. The analysis 
of the AF-HARQ protocol is also quite close to the simulation results. The small deviations 
in the AF-HARQ case is due to the pdf approximations in (|34l) and in (|35] ). Such deviation 
increases as the relay terminal deviates from the center location. However, it is at most several 
percentages as shown in Figures \5\ and [6l 

1 ) Performance comparison: Figure [5] shows the maximum throughput of the cooperative 
HARQ protocols according to the long-term averaged SNR (Es/Nq) when d = 0.5. From this 
result, we can compare the performance of the cooperative HARQ protocols. One observation is 
that DFl and DF3 outperform DF2 over all ranges of SNR due to the packet combining at the 
destination terminal and the maximum throughput of DF3 is slightly better than that of DFl. 
Among the AF-HARQ protocols, the maximum throughput of AF2 is greater than that of AFl. 
Under the peak power constraint, the average transmission energy in Step 2 of DF3 and AF2 is 
greater than that of DFl and AFl, which respectively results in better performance of DF3 and 
AF2 over DFl and AFl. Interesting observation can be obtained by comparing the DF-HARQ 
protocols with the AF-HARQ protocols. Previously it was believed that the AF protocol has the 
same performance as the DF protocol in terms of the D-M tradeoff [i9j| . However, combined with 
the HARQ schemes, the DF-HARQ protocols outperform the AF-HARQ protocols over all range 
of SNR when d = 0.5. The reason is that, in the AF-HARQ protocols, the relay cannot transmit 
additional redundancy because the relay does not decode but just forwards the amplified packet. 
However, in the DF-HARQ protocols, the relay can transmit additional redundancy. Thus, the 
performance of the DF-HARQ protocols is better than that of the AF-HARQ protocols. 

Figure [6] shows the maximum throughput of the cooperative HARQ protocols according to the 
relay location (d) when the long-term averaged SNR is 4dB and 12dB. Among the DF-HARQ 
protocols, DF3 outperforms DF2 over all ranges of d and outperforms DFl in the region of 
d < 0.5. When d > 0.5, DFl and DF3 show similar performance in which the space-time code 
does not improve the performance because the source signal is much more attenuated than the 
relay signal. Also, it is observed that although DFl outperforms DF2 over all ranges of d when 
the SNR is high (12dB), DFl can be worse than DF2 when the SNR is low and d is small due 
to the high error probability at the first transmission in DFl. Among the AF-HARQ protocols, 
AF2 outperforms AFl over all ranges of d and the performance difference slightly increases as 
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the relay tenninal deviates from the center region. The DF-HARQ protocols also outperform the 
AF-HARQ protocols over all ranges of d. Also, it is observed that the performance of DFl, API, 
or AF2 increases as the relay terminal gets closer to the center location, while that of DF2 or 
DF3 does as the relay terminal gets closer to the source terminal due to the space-time coding. 

2) Delay constraint: Figure |7] shows the maximum throughput of the cooperative HARQ 
protocols according to the relay location (d) when M is 2. Although the maximum throughput 
of each protocol decreases, the relative performance among the DF-HARQ protocols (or the 
AF-HARQ protocols) is similar in the case when the maximum number of transmission is large 
(M = 20). However, in this delay constraint case, all DF-HARQ protocols cannot outperform 
the AF-HARQ protocols, since the performance gain by additional redundancy from the relay 
terminal of DF-HARQ protocols is limited. When the relay is not near the destination terminal, 
DFl and DF3 are better than the AF-HARQ protocols. When the SNR is low (4dB), AF2 
outperforms all DF-HARQ protocols in the region where the relay terminal is near the destination 
terminal. When the SNR is high, the AF-HARQ protocols outperform DF2, except the region 
where the relay terminal is close to the source terminal. The value of d, which maximizes the 
throughput of each protocol, may also change according to the delay constraint. Although DF2, 
AFl, or AF2 shows its best performance when the relay is near the center, the value of d 
maximizing the throughput of DFl or DF3 moves to the source terminal as the delay constraint 
becomes tight. 

B. Step power constraint 

1 ) Performance comparison: Figure [8] depicts the maximum throughput of the cooperative 
HARQ protocols under the step power constraint. It is observed that it's different to that of 
the peak power constraint. Among the DF-HARQ protocols, DFl and DF3 outperform DF2 
over almost all regions of d, but the maximum throughput of DFl is greater than that of DF3 
oppositely to the peak power constraint case. Similarly, the maximum throughput of AFl is 
greater than that of AF2 in the region of d < 0.6. On the other hand, the maximum throughput 
of AF2 is greater than that of AFl in the region of d > 0.6. The above observations are very 
interesting because they are different from the previously known results of comparing cooperative 
protocols without HARQ. Among the AF-HARQ protocols, AFl and AF2 without an HARQ 
scheme are, respectively, matched to LTW-AF and NAF in [9]. The DM-tradeoff of the NAF is 
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also better than that of LTW-AF. However, AFl can outperform AF2 in many cases. The reason 
is that the transmission power is divided to the source and the relay terminals to obtain diversity 
gain in AF2. However, performance loss by the path-loss is greater than the diversity gain in 
many cases. 

2) Delay constraint: Figure |9] shows the maximum throughput of the cooperative HARQ 
protocols according to the relay location {d) when M is 2. The maximum throughput of each 
protocol degrades compared to the case of M = 20. However, similarly to the peak power 
constraint, the relative performance among the DF-HARQ protocols (or AF-HARQ protocols) 
is similar when the maximum number of transmission is large (M = 20). In the step power 
constraint, all the DF-HARQ protocols cannot outperform the AF-HARQ protocols. The AF- 
HARQ protocols outperform DF2 over all ranges of d. When SNR is low (4dB), AF2 outperforms 
all DF-HARQ protocols in the region where the relay terminal is near the destination terminal. 

C. The initial transmission rate 

In this subsection, the tradeoff between the initial transmission rate and the delay, such as 
the maximum number of transmissions and the average delay, is shown for DFl under the peak 
power constraint. Although not shown explicitly, similar results and discussion can be obtained 
for other protocols. Figure [10] shows the average delay of DFl according to the maximum number 
of transmissions for various initial transmission rates (r) when the relay location {d) is 0.5 and 
the long-term averaged SNR is 12dB. The average delay is given by 

E[T] 

Tavg = ■ (38) 

— Pout 

When the initial transmission is given and the maximum number of transmissions (delay con- 
straint) is smaller than a certain threshold, the average delay increases very sharply (i.e. the 
outage probability goes to 1). Otherwise, when the maximum number of transmissions is greater 
than the threshold, the average delay is almost constant regardless of the maximum number 
of transmissions. Thus, in order to lower the average delay, it is useful to reduce the initial 
transmission rate, but it is not necessary to increase the maximum number of transmissions if it 
is greater than a certain value. 

To compare the throughput among the various maximum number of transmissions (M), the 
worst-case coding-rate, R = r/M, is used. Figure [H] shows the throughput of DFl according to 
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R under the peak power constraint. As can be seen in Figure [HI when R is small, the throughput 
increases as R increases. However, when R is large, the throughput decreases as R increases. 
This is because, the outage probability (pout) increases as R increases and it becomes a dominant 
factor of performance degradation when R too is large. This trend of DFl is similar to that of 
IR-HARQ protocol which is shown in [11811 . Additionally, it is seen that the region of R, where 
the throughput decreases sharply, is consistent over all values of M. Thus, for a given M, the 
constraint for the initial transmission rate can be set to 

r < R*{d,SNR)M, (39) 

where R*{d, SNR) is defined as the infimum of the worst-case coding-rate yielding zero through- 
put when M goes to infinity, which is a function of the relay location d and the long-term 
averaged SNR. Also, if M is sufficiently large, we can reduce the average delay without 
sacrificing the throughput as can be seen in Figure \T0\ and [TTJ However, if M is small, there is 
a trade-off between the average delay and the throughput according to the initial transmission 
rate. 

Figure [121 shows the optimum initial transmission rate in terms of the throughput according 
to the maximum number of transmissions (M) under the peak power constraint. From this 
figure, it is easily seen that the optimum initial transmission rate increases almost linearly as 
the maximum number of transmissions increases. Thus, the optimum initial transmission rate 
for M, r*(M, d, SNR), can be obtained as 

r*(M, d, SNR) = a{d, SNR)M, (40) 

where 

aid, SNR) = mn rW fNR)^ 

In order to apply the above optimization and tradeoff, R*{d, SNR) and a{d,SNR) should be 
easily calculated from the knowledge of d and SNR, which remains as future work. 

VI. Conclusion 

In this paper, three DF-HARQ protocols and two AF-HARQ protocols were shown and ana- 
lyzed using unified frameworks based on the state transition diagram. A unified framework was 
proposed for each type of cooperative HARQ protocols, which can provide actual performance 



December 16, 2008 



DRAFT 



JOURNAL OF MbX CLASS FILES, VOL. 1, NO. 13, MAY 2008 



20 



evaluation with respect to channel and environment. The maximum throughput of each protocol 
was searched under the peak power and the step power constraints. By obtaining real performance 
comparison among protocols, it was shown that the throughput and relative performance of the 
cooperative HARQ protocols varied according to the relay location, the maximum number of 
transmissions (or the delay constraint), the initial transmission rate, and the power constraint. 
Interesting observations from the analysis as follows: i) the maximum achievable throughput of 
the DF-HARQ protocols can be much greater than that of the AF-HARQ protocols due to the 
IR transmission at the relay terminal; ii) protocols having worse D-M tradeoff without HARQ 
scheme can outperform protocols having better D-M tradeoff in many cases, when HARQ scheme 
is combined; and iii) the region maximizing the throughput of each protocol changes according 
to the relay location, the maximum number of transmissions, and the long-term averaged SNR. 

Additionally, the tradeoff between the initial transmission rate and the delay constraint (or 
the average delay constraint) was observed. However, if the maximum number of retransmission 
is sufficiently large, we can reduce the average delay without sacrificing the throughput. Also, 
it was shown that there is an optimum initial transmission in terms of throughput and it is 
almost a linear function of the maximum number of transmissions, in which the gradient is the 
function of the relay location, and the long-term averaged SNR. Developing a simple algorithm 
that evaluates the optimum initial transmission rate for a practical system remains future work. 

Appendix A 
The probabilities for each state are given as 

v4i=X„ (41) 
A^=am-iAm-i, (42) 

Rl,n bn-^m -^m,n 'J^m— l,n-Rm— l,n) (43) 
M M 

■^e ^ ^ Cj72^m ~l~ ^ ^ (im,nRm,n- (44) 
m=l m=l 

By substituting the transition probabilities shown in Table H] into (|42)) - (|44l) . we obtain 

A^=p^ (m - 1) p2 (m - 1) T'"-iAi, (45) 
Ri,n=Pi {n) {p2 (n - 1) - p2 (n)) TMi, (46) 

Rm,n=P3,n - 1) T"^^^i?l,„. (47) 
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Using the above equations, can be simplified as 

/ M 

Xe^Ai I qi (m)p2 (m - 1) T^'U'^X + Z 



M 



+ ^1'- S «3.- {m)T^U' + Z'), (48) 

n=l \"^=1 ^ 

where Z ^ pi (M) p2 (M) and Z' = p3,„ (M) T^. The transfer function of the state transition 
diagram of the DF-HARQ protocols is then given by 



s 

M 



m=l 

+ ( E^^i (^) 52 in) T^'g (T, A)] , (49) 



where 



^n=l 

M 



g{T,A)^[J2 ?3,n (m) T^C/'- + ps,n (M) 



^,m=l 

The average transmission number E [T] and the average successful transmission rate E [TZ] can 
be obtained by differentiating the transfer function by T and U, respectively. Thus the average 
transmission number is 

df (T, A) 

T=1,A=1 



E[Th 



dT 

M 

=^ mqi (m) p2 (m - 1) + Mpi (M) p2 (M) 



m=l 



+^npi (n) 52 Hc/(1, 1) +Pi (n) ^2 ^^'(l, 1), 



n=l 



(50) 



where g'{l, 1) = (^^^iiTiqs^n {''tt) + Mp^^n [M]^. It can be further simplified as 

M-l M 

£; [T] = 1 + ^ pi (m) p2 (m) + J^Pi (n) {n) E [T% (51) 



m=l n=l 
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where E [T'] = 1 + J2m=i P3,n ("^)- From the transfer function in (|49l) . the average successful 
transmission rate is obtained as 

dfiT,A) 



E[n]- 



(52) 

T=1,A=1 



dA 

--r (1 - Pout) , (53) 



where 

M 



Pout = Pi Pi {f^) 12 {m) P3,m {M) (54) 

m=l 

is the outage probability of the DF-HARQ protocols. By substituting (ISTI) and (|52|) in the 
throughput of the DF-HARQ protocols can be obtained. 

Appendix B 

In the state transition diagram of the AF-HARQ protocols, the probabilities of each state are 
given as 

Ai=Xs, (55) 

Am=bmRm.-l^ (56) 

Rrn (^m-^mi (57) 
M 

X(, = ^^ CmAm + dmRm- (58) 
m=l 

By substituting the transition probabilities shown in Table HI] into (|56l) and (|58l) . we obtain 

Am=Pi{m - l)p2(m - 1)T2(™-i)Ai, (59) 
i?^=Pi(m)p2(m - 1)T2— 1^1, (60) 

Using the above results, X^. can be simplified as 

M 
m=l 

+ Pi(m)g2(m)r2'"f/'-X,. (61) 
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Thus the transfer function of the AF-HARQ protocols is given by 

M 

= J2 qi{m)p2{m - l)T^"^-^U' 

m=l 

+Pl(m)g2(m)T2™[/^ (62) 

By differentiating the transfer function by T and U, the average transmission number and the 
average successful transmission rate can be respectively obtained as 



E[T]- 



df{T,U) 



\T=l,U=l 



dT 

M-l 

=1 + ^ pi{m)p2{m) 



m=l 
M 

+ ^Pi{m)p2{m- 1), (63) 

m=l 

and 



df{T,U)^ 



\T=1,U=1 



dU 

--r{l-pi{M)p2{M)). 



Thus the outage probability of the AF-HARQ protocols is obtained as 

Pout=Pi{M)p2{M). (64) 
By substituting (|63l) and (|64|) in (O, the throughput of the AF-HARQ protocols can be obtained. 
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Fig. L The state diagrams of the cooperative DF-HARQ protocols. 
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Fig. 2. The state diagrams of the cooperative AF-HARQ protocols. 
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Fig. 3. The state transition diagram for the DF-HARQ protocols. 
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Fig. 4. The state transition diagram of the AF-HARQ protocols. 
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Fig. 5. Maximum throughput comparison among protocols according to SNR under the peak power constraint. 
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Fig. 6. Maximum thiroughput comparison among protocols according to d under the peak power constraint when M — 20. 
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Fig. 7. Maximum tliroughput comparison among protocols according to d under the peak power constraint when M — 2. 



M=20, step power constraint 
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Fig. 8. Maximum throughput comparison among protocols according to d under the step power constraint when M = 20. 
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M=2, step power constraint 
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Fig. 9. Maximum throughput comparison among protocols according to d under the step power constraint when M = 2. 



DF1, Es/No=12dB, d=0.5 



20 



t 15 



i 10 

> 
< 





r=10 




r=20 




r=50 




r=100 



20 40 60 80 

The maximum number of transmissions (M) 



100 



Fig. 10. The average delay of DFl according to the maximum number of transmissions (M) under the peak power constraint. 
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Fig. IL Throughput of DFl according to worst-case coding-rate (R) under the peak power constraint. 




Fig. 12. The optimum initial transmission rate according to the maximum number of transmissions (M) under the peak power 
constraint. 
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